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Abstract 



We consider the potential problems due to the production of inflatinos and gravitinos after in- 
^ ■ nation. Inflationary models with a single scale set by the microwave background anisotropics 

have a low enough reheat temperature to avoid problems with the thermal production of 
gravitinos. Moreover, the nonthermal production of gravitinos has been shown to be suffi- 
ciently small if the sector ultimately responsible for supersymmetry breaking is coupled only 
gravitationally to the inflationary sector. Still, in some models, inflatinos can be created 
during preheating with a substantial abundance. The main contribution to the gravitino 
abundance may thus come from their decay into the inflaton, or into its scalar partner, as 
S^ well as from the inverse processes. We show that this production needs to be strongly sup- 

pressed. This suppression can be realized in the simplest scenarios which typically have a 
sufficiently high inflationary scale. 
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1 Introduction 

It is well known that in the particle spectrum of the minimal supersymmetric standard model 
(the MSSM with R-parity conservation), the lightest stable particle (LSP) is an excellent 
cosmological dark matter candidate fl], [|. Indeed, in the constrained version of the MSSM 
based on minimal supergravity, one can easily satisfy all of the recent accelerator constraints, 
and at the same time have an LSP with an acceptable and significant cosmological relic 
density ||. However, it is also well known that in these same models, there is a potential 
problem with the abundance of gravitinos in the early Universe Q. 

Inflationary models are also subject to potential gravitino problems. If the Universe re- 
heats to sufficiently high temperatures, gravitinos will be produced thermally with excessive 
abundances || ||, || . For a gravitational decay of the inflaton field the reheat temperature 
is directly related to the inflationary scale. When this scale is fixed to match the observed 
size of density fluctuations, the reheat temperature is typically low enough to avoid thermal 
overproduction. Recently, the non-thermal production of gravitinos has also been consid- 
ered 0. It has been shown || that in this case too, gravitinos are produced with a sufficiently 
low abundance, provided the inflationary sector and the one responsible for supersymmetry 
breaking are distinct and only gravitationally coupled. 

The nonthermal production is much more relevant for the fermionic partner of the infla- 
ton, the inflatino, which in some cases can be produced with a significant abundance. Thus, 
gravitinos may be overproduced through inflatino decay, if the channel inflatino — > inflaton 
(or its scalar partner) + gravitino is kinematically allowed. Depending on the relative masses 
of the inflaton and inflatino, significant production may be expected instead by the inverse 
process. This is particularly true if the inflaton sector is only gravitationally coupled to 
matter, since in this case the above decays will have a rate comparable to the one generating 
the thermal bath. Indeed, in such a scenario, the decay channels into gravitinos need to 
be strongly suppressed. We discuss the possible kinematic suppression of this channel and 
derive a strong upper bound on the scale of inflation. We further show that this bound is 
satisfied by the simple single scale (supergravity) inflationary models. 

We begin by reviewing the gravitino problem with respect to both thermal and non- 
thermal production. We focus on simple models of inflation with a single scale set by 
the microwave background anisotropics. In section 3, we formulate the inflatino problem 
and discuss solutions for both the gravitational and non-gravitational decay of the inflaton. 
This analysis is generalized in section 4, where we discuss the merits of lowering the scale of 
inflation. In section 5, we consider two specific models of inflation: 1) a model of new inflation 



based on supergravity, and 2) chaotic inflationary models. A discussion and concluding 
remarks are made in section 6. 

2 The Gravitino Problem 

Unless very massive, 7723/2 > 20 TeV, gravitinos would disrupt the successful predictions of 
big bang nucleosynthesis. The argument is relatively simple: because of their late gravita- 
tional decays, gravitinos dominate the energy density of the Universe which becomes matter 
dominated with a Hubble expansion rate given by H ~ m 3 ( 2 T 3//2 /Mp, where Mp is the (re- 
duced) Planck mass. Gravitino decay occurs when the decay rate r 3 / 2 ~ m 3 , 2 /Mp ~ H or at 
To ~ m 3 , 2 /Mp . Subsequently, the Universe is reheated to a temperature Tr ~ p(To) 1//4 ~ 
m 3 , 2 /M P . The limit on the gravitino mass is obtained by requiring Tr > 1 MeV, thus 
allowing big bang nucleosynthesis to occur in a radiation dominated Universe. 

However, since the gravitino abundance is diluted during inflation 0, this simple con- 
straint is significantly altered. In fact, the constraint should be expressed as a function of 
the gravitino mass and abundance after inflationary reheating. The most restrictive bound 
on their number density comes form the photo-destruction of the light elements produced 
during nucleosynthesis JT0| 



n V2 /s < 10-"(100GeV/m 3/2 ) (1) 

for lifetimes > 10 4 sec. The thermal production of gravitinos regenerated after inflation has 
been estimated |2|, ||, ||, [lCJ 



n 3/2 /s ~ (r/7T)(T 3/2 /T 7 ) 3 ~ aN(T R )(T R /M P )(T 3/2 /Tj 3 (2) 

where T ~ q:N(Tr) (Tp/Mp) is the production rate of gravitinos, N is the number of degrees of 
freedom, and the ratio (T3/2/T 7 ) 3 accounts for the dilution of gravitinos by the annihilations 
of particles between Tr and nucleosynthesis. From the gravitino regeneration rate one can 
derive a bound on Tr 

T r < 10 9 GeV(100 GeV/m 3 / 2 ) (3) 

For recent discussions on this bound, see [|TT] . 



The significance of the above constraint depends on the model of inflation. We assume 
for simplicity that inflation is governed by a single scale such that the inflaton potential is 
of the form V(<p) = A 4 P(0/Mp), where P(cj)/Mp) is some suitable function which generates 
inflation. The scale A can be related to the size of large scale fluctuations measured by 



COBE ||12| . In general this relation is dependent in the choice of P , so that different values 
for A can be considered. However, several of the proposed models typically give [H| 



---^-0(100)-^ (4) 

p 1Ott 3 /20 y ' M p 2 w 



which, once related to the COBE results, implies |l ! 

A 2 



M 2 P 



few x 1(T 8 (5) 



Fixing (A 2 /Mp) has several general consequences for inflation |H|. For example, the 
Hubble parameter during inflation becomes H 2 ~ (A A /3Mp), leading to H ~ lCT 7 Mp. The 
duration of inflation is r ~ Mp/A 4 , and the number of e-foldings of expansion is Ht ~ 
(Mp/A 2 ) ~ 10 9 . If the inflaton decay rate goes as V ~ m^/Mp ~ A 6 /Mp, the Universe 
recovers at a temperature T R ~ (YMp) 1 ' 2 ~ A 3 /M P ~ 10 _11 M P ~ 10 8 GeV. This low 
reheating temperature appears to be safe with regards to the gravitino limit @ discussed 
above. 

Recently, considerable attention has been focused on the possible non-thermal produc- 
tion of gravitinos after inflation [|j], ||. Fermionic quanta can be created fl6j in significant 



amounts during the first stages of reheating (preheating). When applied to supergravity 
models, preheating will lead to the non-perturbative production of the fermionic partner 
of the inflaton, the inflatino, and in general, of any other fermion which is strongly cou- 
pled to the inflaton field. If, for example, there is substantial mixing between the inflatino 
and the longitudinal component of the gravitino, the goldstino, preheating may result in an 
overproduction of gravitinos. As is well known, the goldstino is a linear combination of the 
fermionic partners of the scalars responsible for supersymmetry breaking. During inflation 
and the beginning of reheating, supersymmetry is mainly broken by the inflaton implying a 
strong correspondence between the inflatino and goldstino at this early stage. However, this 
correspondence does not necessarily hold at late times, since supersymmetry may be broken 
by other fields in the true vacuum of the theory. In this case, the final gravitino abundance 
will be much smaller than the inflatino abundance ||. Indeed, this is the case in most 
models as it is natural to distinguish between inflation and supersymmetry breaking due to 
the very different energy scales associated with the two phenomena. The relic abundance of 
gravitinos will thus ultimately be related to the strength of the coupling between these two 
sectors in a given model. 

If the scale of inflation, A, is much larger than the other scales in the theory, the nonther- 
mal production of inflatinos can be accurately computed by just considering the inflationary 



sector. Typically, particle creation occurs during the very first few oscillations of the infla- 
ton on a timescale which is the inverse of the inflaton mass, m^ ~ A 2 /Mp . Most of the 
inflatinos are produced with momenta k < m^ , and hence their initial number density is 
approximately given by M 

nf~0.01mj~0.0lf-^j (6) 

The final inflatino abundance is Yg = ng/s where s ~ p 3//4 is the entropy density at the time 
of inflaton decay. For now, let us generalize our previous consideration of inflaton decay, and 
assume only that the decay is given by V (rather than assuming specifically a gravitational 
decay). In this case, if inflaton oscillations begin when the scale factor is R = Rfp, the energy 
density in oscillations is p^ ~ m^Mp(R ( j } /R) 3 and H ~ m^R^/R) 3 ^ 2 . Decays occur when 
H ~ T, or when R = R^ = (m^/T) 2 / 3 .??^. During inflaton oscillations, the inflatino number 
density also scales as R~ 3 so that at the time of decay, rig ~ 0.01 m^R^/Rdtf,) 3 ~ O.Olm^r 2 ~ 
0.0lA 2 r 2 /Mp. The entropy density is also easily computed at decay, s ~ (T 2 Mp) 3 / 4 ~ T|. 
Thus for a massive inflaton, one finds 

— fi?)~- (£)'(£) 

We see, therefore, that a late inflaton decay (low Tr) results in a small inflatino abundance. 
This is due to the fact that the energy density stored in the coherent oscillations of a massive 
inflaton redshifts as matter (p ~ R~ 3 ), so that the quantity n^/p 3 ^ decreases with time. 

For the case of a gravitational decay, if we take A from eq. (j|) and its corresponding 
reheat temperature Tr ~ 10 8 GeV, we find a negligible inflatino abundance Yg < 10~ 19 . For a 
massless inflaton the final abundance will be significantly larger, due to the stronger decrease 
of the inflaton energy density. For example, for a V ~ 4 potential, the energy density of <fi 
redshifts as radiation, and the last factor Tr/A is absent in the final abundance (|7p. In this 
case, the final inflatino abundance may be as large as Yg ~ 10~ 12 . 

The calculation of the nonthermal production of the longitudinal gravitino component 
in models with more than just the inflationary sector || is more complicated. The mixed 
inflatino-gravitino system is very involved, and analytical solutions for the gravitino abun- 
dance are still lacking. A consistent quantization of the system, with an accurate definition 
of the occupation numbers of the fermionic eigenstates, is however available. For the case 
in which the sectors responsible for inflation and supersymmetry breaking are coupled only 
gravitationally, numerical computations show that gravitino production is restricted to a 
safe level. The longitudinal gravitinos are created on a physical timescale of the order m^i , 
with a physical typical momentum comparable with m 3 / 2 ^C A. 



3 The Inflatino problem 

In the previous section, we have shown that if we restrict the inflationary model so that it 
is coupled only gravitationally to the other sectors of the theory, both the thermal and the 
nonthermal production of gravitinos are reduced to a safe level. We have also seen that the 
nonthermal production of inflatinos is model dependent and can be substantial. If this is 
the case, and if the inflatino is heavier than the inflaton, the subsequent gravitational decays 
of the inflatino (to an inflaton and gravitino) could lead to an overproduction of gravitinos 



(other decays of the inflatino in some specific models are discussed in fl7|l). 

On the contrary, when the inflaton is heavier than the inflatino, then the decay channel 
— > cf) + G may be problematic if kinematically allowed || . The density of gravitinos pro- 
duced by inflaton decays is easily estimated. By comparing the number density of inflatons 
just prior to decay, n^ ~ m^Mpi^R^/ R d ^) 3 ~ m^/Mj, to the entropy density just after decay 
s ~ ml jMp , one finds n^/s ~ (m^/Mp) 1 ^ 2 ~ A/Mp. (It is also convenient to write 
n^/ ' s ~ (TMp) 1 / 2 jm$ ~ Tp/m^). If the branching fraction for inflaton decays to <p + G is 
1/N, then the resulting gravitino abundance is 

y = I*£ „ _A_ „ 10-6 (8) 

3/2 s NM P V ; 

in clear violation of the bound ([I]). This is what we call the inflatino problem. We note 
that, for a purely gravitational inflaton decay, the limit from direct production is stronger 
than that due to thermal production. Indeed, by requiring Y" 3/ / 2 < 1CT 13 , eq. (|[) yields the 
constraint A < 10~ n Mp on the inflationary scale, which is stronger than the bound from the 
thermal production of gravitinos. Requiring T R < 10 9 GeV, and noting that for gravitational 
decays, T R ~ (TMp) 1 / 2 ~ A 3 /Mp } one finds A < 10 15 GeV. 

One may hope to relax this difficulty by allowing for non-gravitational couplings of the 
inflaton to matter. For example, if we suppose that the decay rate to matter is given 
by T m ~ g 2 fn ( f ) , then the gravitino abundance will be suppressed, Y 3 / 2 ~ (rG/T m )n^/s 
where To ~ m^/Mp is the gravitational decay rate. In this case, decays occur earlier and 
n 4>l s ~ giMp/m^) 1 ^ 2 ~ gMp/A. The gravitino abundance produced by inflaton decays is 
now 

Y * ~ WmI (9) 

Therefore, the constraint on A from the gravitino abundance becomes A < fewx 10 _4 (y' 1//3 Mp. 
However, for stronger than gravitational decays (g > A 2 /Mp), the constraint from the ther- 
mal production of gravitinos becomes more significant. From Tp ~ gA , one finds the limit 



A < 10 9 GeV/g . The combination of the two bounds from the thermal and the direct grav- 
itino production is weakest when g ~ 10 -15 / 4 . For this value, the limit A < xl0~ 5 Mp is 
found. Thus the simple generic scenario (with a single scale set by the CMB) is not possible 
if the decay of the inflaton to inflatino + gravitino occurs unimpeded. 

Of course, as noted in 0, if it should happen that the decay is kinematically forbidden 
(jm^ — rrti\ < 77^3/2), then the simple single scale inflationary model works very well. Here, 
we note that if the scale of supersymmetry breaking, /1, is significantly below that of infla- 
tion, i.e., |i« A, then even though the decay <fi — > (f> + G is allowed, it will be naturally 
kinematically suppressed with respect to the other decay channels of the inflaton field. This 
will open an allowed window for A even in the simplest scenarios. 

Independent of the details of the decay, the rate will always carry a final state momentum 
suppression factor. The overall decay rate can be written as T ~ {l/m ( j ) )\M.\ 2 {j)/m ( f > ) ) where 
\M\ is the amplitude for decay and the final state momentum suppression factor is 

1/2 

m 2 . ml / A 2 



2(m| + m§ /a ) {ml-m% 2 f \ f 

2p/m rj> = J 1 v — 2 + — V —— A ] ~ — (10) 



recalling that m$ ~ m? ~ A 2 /Mp and m 3 / 2 ~ fi 2 /Mp. Thus in models in which /iCA, 
there will be a significant suppression in the production of gravitinos by either inflaton or 
inflatino decay (note that additional suppression may come from the specific form of the 
amplitude M. as well). 

If we take into account the suppression factor ([H]) in the bound for the direct gravitino 
production by inflaton decay, and we combine it with the limit coming from the thermal 
production, we find that the scale A must lay within the interval 

10l3 iV7tr ~ A ~ 1Ql5GeV ( n ) 

For n ~ 10~ 8 Mp, the lower bound in (|TT]) is about 10 13 GeV. For smaller values of A 
the kinematical suppression factor ( |l0|) is no longer capable of maintaining a low gravitino 
abundance. We notice that the scale A ~ 1CT 4 Mp as given in eq. (||) is within the allowed 
interval. Indeed, including the suppression factor fllTf ) to eq. (H), the gravitino abundance 
produced by decays is now 

v A ^ 2 m3 / 2 in- 14 (io\ 

Fs/2 ~ nm; A*" ~ tva ~ 10 (12) 

and clearly satisfies the bound ([!]). The simplest models of single scale inflation with purely 
gravitational decays and a scale given by (||), therefore, do not suffer from a gravitino prob- 
lem. 

6 



For completeness, we note that when non-gravitational decays are also allowed, the kine- 
matic suppression allows a wide range of values for A (and g) . While the constraint from the 
thermal production of gravitinos is unchanged (gA < 10 9 GeV), the constraint from decays 
becomes A/j, 2 /gNM P < 1(T 13 which for /i = 10" 8 M P becomes A < I0 23 g GeV. 

Before we move on to specific examples of inflationary models, we note that there are 
at least two other potential sources of gravitino production. First, there is the possibility 
for an additional contribution to the direct production of gravitinos by the decay — > 
GG . During inflaton oscillations, the rate for such decays has been estimated to be r 3 / 2 ~ 
A 10 /M 9 |H|. This implies that the abundance of gravitinos will be y 3 / 2 ~ (T 3 /2/T)n^/s ~ 
(A/Mp) 5 . For all of the models being discussed, this is sufficiently small. Second, the 
remaining scalar degree of freedom in the inflaton supermultiplet may also be problematic 
in certain cases. This field, the sinflaton 0', is normally described by the complex direction 
of the full inflaton potential. The complex direction is stable in most models, and thus 
classically it is not excited. However, if lighter than the inflatino, inflatino decay to the 
sinflaton + gravitino may yield another source of gravitino production. The above arguments 
would therefore also apply to these decays as well. In addition, it is possible that quantum 
fluctuations will excite the sinflaton during inflation. Even though the sinflaton is far from 
being massless, quantum fluctuations will lead to <// ~ H |Tj|. However, in generic models 
since the amplitude of sinflaton oscillations is much smaller than the corresponding amplitude 
for inflaton oscillations, their decay into inflatinos is expected to be suppressed by a factor 
of {H/Mpf ~ (A/Mp) 4 . 

4 Lowering the scales of the inflationary sector 

As we have seen, the potential of the inflaton field during inflation is constrained by the 
magnitude of density fluctuations measured by COBE. In the slow roll regime (3 H cf) ~ — V) 
eq. (U) gives 

V^ /4 ~ 0.027 e 1 ^ M P (13) 

where e 2 = M 2 (V'/V) /2 <C 1 is one of the two slow-roll parameters (prime denoting 
derivative with respect to 0) and the suffix H reminds us that the two quantities have to be 
evaluated when the scales measured by COBE left the horizon, about 60 e-foldings before 
the end of inflation. As we have said, in the simplest models of single scale inflation, this 
relation fixes the scale of the inflaton potential A = V H to be about 10 _4 Mp . However, 
models with a much smaller scale and acceptable density fluctuations can be derived. As 
follows from eq. (13), in models with one single field this can be done at the expense of 
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a small e parameter, that is by taking a very flat potential during inflation. Such a flat 
potential may arise more naturally if more scalar fields are present, as for example in hybrid 



inflationary models p0 |. Due to this freedom, most of the above results of the previous 
section have been given for an arbitrary inflationary scale. 

In this section, A 4 denotes the value of inflaton potential at the end of inflation, when the 
reheating stage begins. Due to the slow motion of during inflation, this scale is typically 
very close to Vh • During reheating, the inflaton field oscillates about the minimum O of V, 
with V (4>o) = . As is typical for a massive inflaton, we assume that the quadratic term 
dominates the Taylor expansion of V around 0o • We denote by F the amplitude of the 
inflaton oscillations at the initial time t ~ H~ l ~ Mp/A 2 . While in the above discussions 
the natural assumption F = Mp (originally dubbed primordial inflation [21]) was made, in 
the following we discuss how the results of the previous section are affected when the value 
of F is lower. 

This analysis is particularly simplified by noticing that most of the bounds previously 
discussed are directly related to the inflaton mass m^ rather than to the scale A . For a 
quadratic potential and generic values of A and F one has m^ ~ A 2 /F . 

Let us first consider the case in which the inflaton decays only gravitationally. For a 
quadratic potential, the reheating temperature is given by T R ~ (T M P ) ' ~ Jm^/Mp , so 
that the thermal bound (0) simply gives 

m < 10 12 GeV (14) 

Following the same line of arguments of the previous sections, the inflaton "abundance" 
at the decay time can be estimated to be 



n <t> ^ T R ^ [™* (15) 



p f 777^ V M P 

The abundance of inflatinos produced nonthermally by the inflaton oscillations is instead 
(also evaluated at the inflaton decay) 

%_ „ 2 rrnTp ^ ,m±\W 

pf~ Ml " iU \Mp) Ubj 

We see that the inflaton abundance is always higher than the inflatino abundance. Also 
the abundance of gravitinos produced nonthermally becomes smaller as m^ decreases. In- 
deed, gravitinos will still be mainly produced at the time t ~ rn^} 2 with a typical momentum 
k ~ 7773/2 5 independent of the value of m^ . However, a lighter inflaton implies a longer life- 
time t^ . Since the quantity 713/2//?/ decreases in the time interval 7773/2 < t < r^ , lowering 
777^ will thus decrease the final nonthermal gravitino abundance. 



Direct gravitino production still mainly occurs through the — > <fi G decay, if kinemati- 
cally allowed. This production is reduced to a safe level as long as 

V2^ 1Ql3 ^ „« 

where the kinematical suppression factor m.3/2 /m<f> has been included. 

Rewriting eqs. ( |I^ ) and ( |T7D in terms of A and F, the allowed window for A becomes 



loI3 i4\/s; SAS10l5GeV \/s; < 18 > 



For F = Mp , this coincides with the result ( |Tl|) obtained in the previous section. For smaller 
F , smaller values of A must be considered, as it is obvious from the scaling m^ ~ A 2 /F . 

As before, additional freedom in the choice of parameters is allowed if the inflaton has 
non-gravitational decays. For T ~ g 2 fn ( f ) and for 7713/2 — 100 GeV, the bounds on the 
reheating temperature and on the direct gravitino production become 

m<j> < GeV/g 2 , m4> < 10 28 GeV g 1 (19) 

Also in this case, the bounds can be written just in terms of the inflaton mass. Notice that 
when the inflaton has nongravitational decays, the limit from direct gravitino production 
gives a lower bound (rather than an upper one) on A . This follows from the different 
scalings of the gravitational and the non gravitational rates with m^ . 

5 Models 

Given the above arguments for the post inflationary production of gravitinos through thermal 
and non-thermal effects, we will now apply them to some specific examples. First we consider 
a simple model of inflation based on N=l supergravity [El, |23|. We subsequently consider 



gravitino production in models of chaotic inflation ||24|| . In both examples, we will assume a 
superpotential of the form 

W T = W ($) + /i 2 (S + 13) 
W{$) = A 2 ^(<l>) (20) 

where we set Mp = 1 for convenience. The first term specifies the dynamics of the inflaton 
field (i.e. the scalar component of $) during inflation and the first stages of reheating. The 
second term is known as the the Polonyi superpotential p5| , and is the simplest example of 
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gravitationally mediated supersymmetry breaking involving only one chiral multiplet S . We 
also assume a minimal Kahler potential G = K + In IWtI 2 , with K = $' $ + S* S. In the 
minimum of the theory, the Polonyi scalar takes the vacuum expectation value, (s) = \/3— 1 , 
up to possible fi 2 = /x 2 /A 2 <C 1 corrections which may arise from the interaction with the 
field ■ The parameter (3 is fine-tuned to 2 — \/3 (again up to possible ft 2 corrections) to 
cancel the cosmological constant. 

In the following examples, we will assume that in the limit of \x — > 0, supersymmetry 
is unbroken, that is, supersymmetry is mainly (in some cases completely) broken by the 
Polonyi field. Therefore, we have F^ < (i 2 F s , where, for a minimal Kahler potential, Fi = 
e G l 2 G l = e K > 2 (dW T /d<pi + (p*W T ) and G { = dG/d^. We also have W (<f>) < W (s) p? . 
From the specific form of the Polonyi superpotential one finds F s = V3fi 2 , and the gravitino 
mass will be 

m 3 / 2 ~e K /V (21) 

These relations follow from the fact that in these models, W (0) and F^/A 2 are at most of 
order /i 4 in the global minimum of the theory. 

With these assumptions, the real and the imaginary components of the Polonyi field 

s = (s\ + is 2 ) /a/2 have masses 



m. 



~e*/VV2V3 , m S2 ~e^V 2 V4-2v / 3 (22) 



Analogously, for the inflaton sector we decompose = (0i + «0 2 ) /V% and we denote the 
inflatino field by <p ■ Then, a simple calculation gives (assuming W < fi 2 W" ^ where 
prime denotes derivative with respect to ) 



m <t>i,4>2 = m 0± m 3/2 U-— j +0 (/i' 



m ~4> _ p K /2 

A 2 " 



r-rffl+r**- 



O A 4 (23) 



where Kq and K% are the first two terms in the expansion of the function K = Kq-\-(x 2 K\ + . . . . 
Note that the masses of the inflatino and of the two inflaton components are nearly 
degenerate, their difference being related to the small gravitino mass m 3 / 2 . This is a con- 
sequence of having taken the supersymmetry breaking scale \x much smaller then the infla- 
tionary scale A . We notice that the specific choice of the Polonyi superpotential (p0|) gives 
rn (h ~ m 6 ~ m 4>i > but that the mass differences are always smaller than the gravitino mass 
m 3 / 2 . Hence, the two potentially dangerous decay channels — > + G and — >• + G are 
kinematically forbidden in these models. 
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It is important to stress that the above relations ( p3[ ) are unaffected by the oscillations 
of the Polonyi field about its minimum. Indeed, for 777.3/2 — 100 GeV we have 



■*■ grav 

rn.3/2 



10" 



A 



10- A M P 



(24) 



where r grav ~ 



rrA/Mp is the rate at which the decay <fi —* <fi + G would occur if not 
kinematically suppressed or forbidden. The oscillations of the Polonyi field start at the time 
2 and are eventually damped by the expansion of the Universe. For t > rru/ 2 , their 



m 



3/2 



3/2 



amplitude scales inversely with time. Since the initial amplitude of the oscillations and the 
value of the Polonyi field in the minimum are of the same order (Mp), eq. ( |2"4| ) shows that 
these oscillations can be soon neglected in the calculation of the mass spectrum. 

Actually, in the following we will ignore the evolution of the Polonyi field. This evolution, 
however, leads to a well known problem in cosmology which is compounded in more modern 



theories collectively called the moduli problem. 26 



5.1 New inflation 

A workable example of inflation is provided by 

A 2 



W($) 



2M f 



(<&- 



(25) 



which is clearly of the form advocated in section 2. 

During inflation and the first stages of reheating, the Polonyi scalar is frozen at s = . 
Moreover, since /j<A, the Polonyi sector is completely negligible at this stage. Thus, for 
the moment, we can set \x = and decompose the inflaton field into its real plus imaginary 
components, <p = (0i + ^2) /v2 ■ Inflation occurs while <pi ~ . Near the origin, the real 
direction is relatively flat, the imaginary component has large positive mass, and it is rapidly 
driven to zero. At the classical level, we can set 2 = , since V" is always positive in the 
imaginary direction. Doing so, one finds 

/ \ 2 



4 



n 



V2- 



>i 



4M P 



>0 



(26) 



V2Mp J [ M P Ml y/2Mp 

where 0o = Mp has been set to have a vanishing cosmological constant in the minimum at 
<pi = \plMp . Notice that in the vicinity of the origin, both the linear and quadratic terms 
in 0i cancel when the exponential prefactor is expanded. Inflation in this model occurs 
for 0i < Mp/ (6 v^J • The inflationary scale A which matches the microwave background 
anisotropies is about A/M P ~ 2.6 x 10 -4 , and the spectral index is found to be n s ~ 0.93 . 
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Let us now consider the late behavior of the fields and reintroduce the dimensionless 
quantity jx = /i/A <C 1 . The presence of a nonvanishing Polonyi field slightly modifies the 
potential for and the position of the minimum. Also, the nonvanishing vev of forces 
a small modification to the Polonyi potential. These effects, which are easily computed in 
an expansion series in fi 2 , must be taken into account in the computation of the spectrum 
in the minimum of the theory. There is some freedom in the choice of the parameters, 
since different combinations of 0o an d (3 give a zero cosmological constant. Up to order fjp 
corrections, one finds 

l + ftf + gfi* 

(2-v / 3)-A 2 /2+(^-^-/)/i 4 (27) 



M, 



p 



Mr 



with / and g arbitrary real numbers. 

Decomposing s = (s\ + is2) /v2 , one finds S2 = , while the real components have their 
minima at [] 



) 



l + {f-l)^+(g-f+l-y/3)jM^ 



y/2M P v ' \ * J 2 

<*> = VS-l + ^ + (/ + ^-?W (28) 



V2M J 



again up to order /i 6 corrections. 

In the minimum, the F— terms are given by 



Fs ~ R -2 , v3 . 4 F c 

^-V3^ +— a* , £ 



:! ^^ / i 2 + ^ / i 4 , -f^-v/3/i 4 (29) 



Hence, the inflaton sector provides a negligible (although non-vanishing) contribution to 
supersymmetry breaking in the minimum of the theory. At late times the goldstino can be 
identified with the Polonyi fermion, up to small fi 2 corrections. The masses of the gravitino 



and of the two components of the Polonyi field are given by the general expressions (21) and 



(E2|), so that we have 



m 3/2 



^ 2 - 8 - 1(rlo (lo^v) ^ 



For these minima one has 



s-°« -^= ^ .*-°<n 
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The general relations (|23| ) for the inflationary sector are also respected. More explicitly, 
one finds 



m 4> 

m 4>2 



M P 

A^ 

M P 

A^ 

~M~ P 



J-n/3 



J-n/3 



1+ (/ + f-|)- 



-V3 



!-(/- 2,)^ 



7 



/-- + V / 3)/r- 



(31) 



G and 



+ G are indeed kinematically 



confirming that the two channels 0i — > < 
forbidden. 

With regard to the nonthermal production of inflatinos and gravitinos, we have verified 
by numerical calculations that the model fl2o]) gives results in very good agreement with the 
ones obtained in || for the chaotic inflationary superpotential W ($) = m 2 , $ 2 /2 (see below). 
Inflatinos are produced during the first oscillations of the inflaton field (t ~ m^ 1 = Mp/A 2 ). 
Their final abundance Yi is given in equation (0). As we have remarked, the final value of 
Yi is typically very small, due to the dilution of ni/p^ during the inflaton oscillations in 
the massive inflaton case. The nonthermal production of gravitinos is instead mainly due to 
the Polonyi oscillations, and indeed it occurs on a timescale t ~ m" 1 ~ ^-3/2 • As we have 
already mentioned, their final abundance is much smaller than that of inflatinos, since they 
are typically produced with a physical momentum k ~ m 3 / 2 • 



5.2 Chaotic inflation 

For completeness, we consider two types of chaotic models with superpotentials 

1 X . 

W{$) = - m^ $ 2 and W ($) 



2 .-, 3 -* 3 (32) 

which for <^C Mp lead to the usual quadratic and quartic potentials typical of chaotic 
inflation |M]. As it is well known, the above superpotentials do not lead to inflation in 
supergravity, since corrections from the Kahler potential spoil the flatness of V (0) for > 
Mp (at least in minimal supersgravity). One can still assume that some corrections, relevant 
at high <p , will generate a sufficiently flat potential so to render inflation possible, and that the 
superpotentials fl32l ) will be accurate enough during the reheating stage (see [^7] for example, 
for chaotic models fully consistent with supergravity). After inflation the production of 
gravitinos can be very well estimated with the superpotential (^), at <fi <^C Mp ||. However, 
the two parameters m$ and A cannot be directly linked to the size of the density fluctuations 
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as in eq. (J|), and the "standard" values m^ ~ 10 13 GeV and A ~ 1CT 13 should be taken 
only as indicative ones. In both cases, <fi = is the minimum of the model, and the inflaton 
sector does not contribute to supersymmetry breaking (F^ = 0) at late times. 

Let us first consider the superpotential W = m^ $ 2 /2 . Evaluation of eqs. ( |23"| ) for the 
masses of the inflaton components simply gives 



m i = e2 m <t> > ^i,<^ 2 = e 2 ^ 



, M 2 L ^ 



(33) 



The nonthermal production of inflatinos and gravitinos in this model has been studied in ||. 
The results have been briefly summarized in section Q and at the end of the previous sub- 
section. 

Let us now consider the superpotential W = A $ 3 /3 . In this model as well, the nonper- 
turbative production of inflatinos and gravitinos is analogous to the ones already discussed. 
Inflatinos are mainly produced during the first few oscillations of the inflaton field, with a 
typical momentum k ~ m^ ~ vA0o > 0o ~ Mp being the initial amplitude of the infla- 
ton oscillations. There is however an important difference, due to the fact that, during the 
oscillations, in the V (<p) ~ <f) 4 model the inflaton field has the equation of state of radiation, 
and its energy density scales as R~ 4 . Thus the energy density of <fi oscillations is given by 

P(t> ~ A0* (R+/R)* (34) 

In this phase both ni and s = pj scale as R~ 3 . The dilution factor then drops out from 
their ratio (cf. the last factor in eq. fl7|)) and one is simply left with 

s^M.^af (35 ) 

a (\<f> A f \io- 13 J 

As for the previous cases, nonthermal gravitino production is expected to occur mainly 
during the oscillations of the Polonyi field, with a negligible final abundance. 

In the minimum of V(<p) the masses of the two inflaton components, as well as that of 
the inflatino, vanish, although the inflaton is still massive during its oscillations. Since the 
amplitude (f> changes very slowly with respect to the oscillations timescale, we can get the 



inflaton mass from the average of Jd 2 V/cP(j) during one oscillation. Up to a numerical factor 
of order one, we can set 

rrty(t)~VA0(t) (36) 

It is interesting to note that in this model the inflaton cannot decay gravitationally, since 
with the above mass one has 

^r^A 3 ^/^ 3 
grav "Mj,- M 2 P \RJ [ ' 
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This should be compared to the Hubble rate 

h ~^^{-r) (38) 

Clearly the gravitational decay rate is always smaller than H. For a non-gravitational decay 
with a rate T = g 2 m§ ~ g 2 \ 1 / 2 4>o(R ( p/R), which should be compared to eq. (|38|). Thus the 
scale factor at decay is (R$/Rd<p) ~ Mpg 2 /(po and the inflaton lifetime is given by 

compared here with the timescale of the Polonyi oscillations. Thus although the inflatino 
abundance in eq. fl35|) is somewhat high, its decay (which is necessarily non-gravitational) 
does not lead to gravitino production. 

6 Discussion 

The simple models discussed in this letter were chosen to highlight the potential problems 
induced by the inflatino. While more general remarks have been made in earlier sections, 
explicit computations have been performed for a single inflaton supermultiplet along with 
the supersymmetry breaking sector described by a Polonyi superfield. As remarked earlier, 
we have ignored the cosmological problems associated directly with the Polonyi scalar's 
evolution. We have also neglected the potentially important evolution of other scalar fields 
or other sources of entropy production. 

Since the MSSM contains many flat directions, the cosmological evolution along these 
directions can have a direct impact on the issues discussed here. For example, it is well 
known that some flat directions may be responsible for generating the baryon asymmetry of 



the Universe 2e . In Affleck-Dine baryogenesis scenarios, if the initial field values of the flat 



direction is larger than a critical value, 770 > m 3 / 2 M P jm! , where r\ is the A-D condensate, 
the dominant source of entropy comes from 77 decays rather than inflaton decays, although 
in this case, the baryon asymmetry is generally too large. For smaller initial values of r] , the 
A-D condensate decays during inflaton oscillations. Depending on the value, 770, late entropy 
production by moduli decay may be a welcome feature to bring the baryon asymmetry down 

We have shown that the decay of inflatons to/from inflatinos + gravitino is potentially 
a serious problem for a generic inflationary model based on supergravity. This problem is 
particularly enhanced when the inflaton sector is coupled only gravitationally to matter. In 
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this case the decay channels with a final state gravitino can have rates comparable to the 
decays into matter. This existence of these gravitino producing decays depends on the mass 
spectrum of the model. In the particular case of Polonyi-like supersymmetry breaking, we 
have shown that such decays are in fact kinematically forbidden, since the mass difference 
between the inflaton and the inflatino is smaller than the gravitino mass. 

This fortunate property is not expected to hold for a general supersymmetry breaking 
potential. More generally, the direct gravitino production must be forbidden by imposing a 
strong hierarchy between the scale of inflation A and the size of supersymmetry breaking. 
In gravitational mediated supersymmetry breaking schemes (7713/2 ~ 100 GeV) this can be 
translated into an upper bound on A . This limit weakens as the initial size of the inflaton 
oscillations is decreased. For initial oscillations of the order the Planck scale (primordial 
inflation) the bound is quite severe, and a safe value for the inflationary scale is A > 
10 13 GeV. Models of single scale inflation with a scale dictated by the size of the density 
fluctuations are typically compatible with this bound. 
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